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Abstract. We present results of VLT/FORS2 spectroscopy of galaxies at z ~ 3 in the Hubble Deep Field - South 
(HDF-S). A sample of galaxies was drawn from the photometric redshift catalogue based on the HST/WFPC2 
optical images and the deep near-infrared images obtained with VLT/ISAAC as a part of the Faint Infrared 
Extragalactic Survey (FIRES) project. We selected galaxies with photometric redshift between 2.5 and 4. Most 
of the selected galaxies are bright in rest-frame UV wavelengths and satisfy color selection criteria of Lyman 
break galaxies (LBGs) at z ~ 3. The number of target galaxies with 7ab < 25.0 was 15. We identified new 5 
firm and 2 probable redshifts in addition to confirmations of previously known 6 galaxies at z ~ 3. We found 6 
among these 13 galaxies lie at a quite narrow redshift range at z = 2.80 ± 0.01. Their spatial distribution is fairly 
concentrated and is at the edge of the HDF-S field, suggesting the possible existence of larger galaxy clustering. 
We examined stellar populations of the galaxies with spectroscopic redshifts through comparisons of their optical 
and near-infrared photometry data with template spectra generated by a population synthesis code. The ages 
from the onset of star formation for these star-forming galaxies with I < 25.0 are typically 50-200 Myr, and 
their stellar masses are between (0.5-5) x10 10 Mq, consistent with previous studies. We also compared these SED 
fitting results with those for "distant red galaxies" (DRGs) at z > 2 discovered by FIRES. DRGs have larger 
stellar masses, larger dust attenuation than our UV-luminous LBG sample, and their star formation rates are 
often comparable to LBGs. These trends suggest that majority of DRGs are indeed the most massive systems at 
the redshift and are still in the active star-forming phase. Unless the number density of DRGs is much smaller 
than LBGs, estimates based on UV selected sample could miss substantial part of stellar mass density at z ~ 3. 

Key words, cosmology: observations - galaxies: high-redshift - galaxies: distances and redshifts - galaxies: evolution 
- galaxies: stellar content 



1. Introduction 

In recent several years, there has been much progress 
in understanding the nature of galaxies at high-redshift. 
Deep imaging of blank fields and the photometric redshift 
technique based on multi-band photometry have enabled 
a construction of large samples of high-redshift galaxies 
which are so faint that we need long exposures to obtain 
good signal-to-noise spectroscopy of them. The Lyman 
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break technique, a simple variant of the photometric red- 
shift techniques, intends to detect high redshift galaxies 
by finding a spectral break at redshifted Lyman limit and 
the discontinuity at shortward of Lya due to attenuation 
by intergalactic neutral hydrogen. The number of galaxies 
selected by the technique, so-called Lyman break galaxies 
(LBGs), amounts to an order of 10 3 at z ~ 3 (e.g., Steidel 
et al. ITOOOl 120051 Foucaud et al. l2005|) . and several hun- 
dreds even at z ~ 5 (Iwata et al. 120051 Ouchi et al. l2004"|) . 
Follow-up studies on LBGs such as optical spectroscopy 
(e.g., Teplitz et al. 120001 Pettini et al. I200T1 Shapley et 
al. 120031 Ando et al. 2002J and analysis of their stellar 
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populations using optical-to-near-infrared spectral energy 
distribution (SED) (e.g., Sawicki and YeeEP98, Shapley 
et al. [2001. Papovich et al. I20U1|1 have been made. By 
these studies basic properties of LBGs have been clarified; 
their star-formation rates amount to several tens to sev- 
eral hundreds Mq /year, dust attenuation is fairly small 
(E(B — V) < 0.5), their stellar populations are dominated 
by relatively young massive stars (ages less than 1 Gyr in 
most cases) and stellar mass is an order of 10 10 Mq (about 
1/10 of the Galaxy). 

Although these studies have clearly indicated that the 
Lyman break technique is a powerful method to probe 
the star-forming galaxies at z J> 3, the number fraction of 
LBGs among the total populations of galaxies at a red- 
shift considered remains unclear. Since the Lyman break 
technique selects galaxies which show prominent discon- 
tinuity at A = 1216A and have flat UV spectra, sample 
galaxies are biased to actively star-forming galaxies with 
relatively small amount of dust attenuation. Recently, sev- 
eral studies based on deep near-infrared imaging revealed 
substantial populations of massive galaxies which are faint 
in optical bands but bright in near-infrared (e.g., Franx et 
al. 12777751 Daddi et al. 1277771 Chen and Marzke l277774) . and it 
has been suggested that optical (i.e., rest-frame UV) based 
samples may miss a significant amount of massive galax- 
ies. However, due to their optical faintness, only a small 
fraction of these red galaxies have spectroscopic redshifts. 
Thus current estimates of their redshift distribution and 
number density have large uncertainties. Moreover, since 
the red color can be explained by both heavy dust at- 
tenuation and old stellar population, the nature of these 
galaxies - dusty star-forming galaxies or massive galaxies 
dominated by old stellar population - is unresolved. 

The Hubble Deep Fields (HDFs) are one of the 
most prominent observations for the study of high-z 
galaxies. Deep imaging covering the optical wavelengths 
by HST/WFPC2 (Williams et al. Casertano et 

al. [2000 ) as we H as near-infrared images obtained by 
HST/NICMOS and ground-based observatories enabled 
the photometric redshift estimation down to Jab = 27.5, 
for both northern and southern HDFs. (e.g., Fernandez- 
Soto et al. IE77777I Yahata et al. 12077711 Labbe et al. 12717781 
(hereafter L03); Vanzella et al. 1277774) . In the HDF-S we 
carried out a deep 011+44 narrow-band imaging using 
VLT/FORS1 to examine the escape fraction of Lyman 
continuum from galaxies at z ~ 3 (Inoue et al. 2005). 
Although precise redshift measurements were required to 
use the narrow-band image to discuss the Lyman contin- 
uum escape fraction, in the HDF-S only a handful spec- 
troscopic observations have been made so far (Cristiani et 
al. 120001 Vanzella et al. 12777721 Sawicki and Mallen-Ornelas 
2003 Rigopoulou et al. 2005 ) , and in particular, the num- 
ber of galaxies identified to be at z > 2.5 is still quite lim- 
ited; only 2 galaxies at z > 3 and 8 galaxies at 2.5 < z < 3 
have been reported. In order to increase the number of 
sample galaxies, we carried out spectroscopic observations 
of galaxies with photometric redshifts estimated to be at 
z ~ 3 with VLT/FORS2. The precise redshift information 



is also useful to reduce the uncertainty of physical prop- 
erties (stellar mass, ages, amount of dust attenuation) in 
SED fitting. In this paper we report the results of the 
spectroscopic observations. In section 2 we describe the 
sample selection, observations and data reduction. The 
results and a discussion on a spatial distribution of sam- 
ple galaxies are presented in section 3. We also investigate 
the stellar population and star formation history of the 
galaxies (section 4) , and compare them with properties of 
"distant red galaxies" discovered by Franx et al. (2003) 
in section 5. A summary and conclusions of the paper are 
presented in section 6. 

Throughout the paper all magnitudes are presented in 
AB system. We adopt a set of cosmological parameters of 

(n Mj n A> h) = (0.3, 0.7, o.7). 

2. Observations and data reduction 

2.1. Sample selection and observations 

The target field is the WFPC2 main field of the Hubble 
Deep Field - South, centered at R.A. = 22:32:56.2 
Decl=-60:33:02.7. The WFPC2 observations, data reduc- 
tion, data quality and object detection procedures are 
described in Casertano et al. I|2000[l . Ultra-deep near- 
infrared Js, H, Ks-hand imaging observations using 
VLT/ISAAC have been carried out as a part of the Faint 
Infrared Extragalactic Survey (FIRES). L03 utilized their 
near-infrared images and publicly available HST/WFPC2 
optical images to estimate photometric redshifts for faint 
galaxies. The catalog contains 134 entries with a photo- 
metric redshift between 2.5 and 3.7. We selected a primary 
sample for spectroscopic observation with the Isu cor- 
rected isophotal magnitude brighter than 25.0. The num- 
ber of galaxies satisfying the criteria was 23. 

Figure shows the distribution of sample galaxies in 
the U300 — -B450 versus -B450 — Igu two-color diagram. For 
those without U300 detection, lower limits of Usoo — -B450 
colors are shown. Since the selection of the sample galax- 
ies were not based on the Lyman break method but on the 
photometric redshift using multi-band imaging data, not 
all the sample galaxies satisfy the criteria for selection of 
z ~ 3 LBGs defined by Madau et al. (1996). All except two 
objects with spectroscopic redshift determinations (shown 
with filled symbols in figure^) satisfy the LBG color crite- 
ria (see §3). Two objects with spectroscopic redshifts lying 
at the outside of the LBG color area are HDFS0001 and 
HDFS0085. These two objects locate close to the edges of 
the WFPC2 image, and photometric errors in [/3oo-band 
are estimated to be ~ 0.5mag. Thus we expect that prop- 
erties of sample galaxies with spectroscopic redshifts are 
quite similar to LBGs. 

We made a multi-object spectroscopy using 
VLT/FORS2 in MXU mode. In this mode exchange- 
able mask unit with slitlets prepared on a focal plane 
were used. We defined the positions of slitlets on the 
HST/WFPC2 7 8 i4-band image, using the software pro- 
vided by ESO. For some galaxies with Iab < 25.0, we 
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Fig. 1. L/300 — -B450 versus -B450 — Igi4 diagram for sam- 
ple galaxies for VLT/F0RS2 spectroscopy. The galax- 
ies for which we successfully identified their rcdshifts are 
represented with filled symbols, and those observed but 
failed in redshift determination are shown as open sym- 
bols. Squares are for galaxies which have been detected in 
J73oo-band, and triangles are for those without [/300-band 
detection. Larger marks and smaller ones are for galax- 
ies with J 814 < 25.0 and those with 25.0 < I 8 i4 < 25.5, 
respectively. We also show the color tracks of model star- 
forming galaxies with and without dust attenuation as 
two solid lines. The model was computed assuming a con- 
stant star-formation rate and an age of 100 Myr. The left 
track is for the model without dust, and the right one is 
with dust attenuation (E(B — V) = 0.4). These lines in- 
dicate how the colors of a star-forming galaxy change as 
its redshift increases. The dashed and dot-dashed curves 
are calculated using template spectra of elliptical and Sbc 
galaxies, respectively, from Coleman et al. |1980|). We put 
open circles at redshift interval of 0.5. Points for z > 2.5 
were enlarged. The color selection criteria used by Madau 
et al. I|1996[) for selecting LBGs at z > 2 are indicated by 
a dashed line. 



could not set a mask due to overlappings with other 
target objects. On the other hand, we set some additional 
slitlets for fainter objects, i.e., 25.0 < Iab < 26.0, if 
they did not conflict with slitlets for bright objects. 
The final sample for spectroscopy consists of 15 galax- 
ies with 22.9 < Iab < 25.0 and 10 galaxies with 
25.0 < Iab <S 26.0. In table ^ basic parameters of the 
sample galaxies are listed. 

The observations were executed in service mode, on 7 
nights during 25th June to 28th July 2003. We used the 
300V grism, which provides a 1.6A / pixel resolution over 
3300A- 6600A. Slit widths were 1.7", and the resolving 
power at 5900 A was 260. Exposure time for each shot was 
999 sec in most cases, and the total on-source exposure 
time was 6.1 hours. 



2.2. Data reduction 

Basic data reduction steps, namely, bias subtraction 
(using overscan regions), cosmic-ray removal and flat- 
fielding have been made using IRAF. After the registra- 
tion of all frames, the images were combined using the 
IRAF/IMCOMBINE task. The sky subtraction and ex- 
traction of one-dimensional object spectra were made us- 
ing an interactive program developed by one of us (I.I.). 
Using that program a user can make a background fit- 
ting in any wavelength range and preview the result of 
subtraction before execution. 

The flux calibration and sensitivity correction were 
made using spectral data of standard stars (Feige 110, 
LTT 1020 and LTT 6248) taken during the observing runs. 

3. Results of spectroscopy 

3.1. New redshift identifications 

We firmly identified 5 new spectroscopic redshifts. All of 
these new redshift determinations are based on the iden- 
tifications of three or more spectroscopic features, such as 
Lya emission, low- ionization absorption lines (e.g., Si 11 
and C 11), and the discontinuity of continua at 1216A. 
They are labeled with asterisks in the sixth column of 
tabic H ln figure and |21 we show the extracted one- 
dimensional spectra of galaxies with redshift identifica- 
tions. There are six galaxies with previously determined 
redshifts; one object HDFS00271 is a faint (I = 25.45) 
Lya emitter, and the others are I < 24.4. Redshifts for 
these objects listed in table are based on our spec- 
troscopic data, and the differences from previous studies 
are less than 0.003. There is also one galaxy HDFS 10652 
(HDFS 223250.05-603356.8) which we did not observe 
but the redshift has been determined to be z = 2.652 
by Vanzella et al. I|2002[) . Thus in total 12 galaxies with 
2.5 < z < 3.5 have been found in the HDF-S field. There 
are also two objects whose redshifts are likely to be newly 
determined by our observation, but the identifications are 
not definite because the S /N of spectra was poor and only 
one or two lines are firmly identifiable. We list these two 
probable redshifts in tableHwith parentheses. For galaxies 
whose spectra have Lya line as emission lines we show the 
rest-frame FWHMs and equivalent widths (EWs) in table 
□ For HDFS00565 and HDFS00526 we measured FWHMs 
and EWs for combined spectra, because their positions 
are only 2."1 apart and spectra are blended. We tried to 
extract spectra of their central parts separately and saw 
no difference in their redshifts (z = 2.789). Because the 
signal-to-noise ratios of spectra are poor, we could not 
measure EWs of absorption features. 

The spectroscopic redshifts obtained for our sample 
galaxies are always smaller than photometric redshifts 
estimated by L03. In Fig. 0] we show a comparison of 
spectroscopic and photometric redshifts for all galaxies in 
the HDF-S which have a published spectroscopic redshift. 
There is no object at z > 2 for which photometric red- 
shift is equal to or smaller than spectroscopic redshift. On 
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Spectroscopic Redshift 

Fig. 4. Comparison of Spectroscopic redshifts (z spcc ) and 
photometric redshifts (z p hot) estimated by L03. The differ- 
ence between the spectroscopic redshifts and photometric 
ones are shown along with the spectroscopic redshift. The 
error bars represent the differences (z spe c — Zphot) in cases 
of using lower and upper limits (68% confidence level) of 
photometric redshifts. The galaxies with firm redshifts in 
our study are shown with filled circles. Open circles rep- 
resent galaxies listed in L03. Galaxies for which z spcc have 
been measured by Sawicki and Mallen-Ornelas ( 2003 ) arc 
shown with filled squares. 



the other hand, the photometric redshifts of galaxies at 
z < 2 are systematically smaller than spectroscopic red- 
shifts. This trend can be also seen in figure 6 of L03 but it 
is not clear due to smaller number of galaxies with spec- 
troscopic redshifts. 
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Fig. 5. Spatial distribution of galaxies with spectroscopi- 
cally confirmed redshifts. Galaxies within a redshift range 
2.789 < z < 2.806 are shown as filled circles. The other 
galaxies are shown as open circles except HDFS00085 at 
z = 3.f6 which is marked as a triangle, fn the spectrum 
of HDFS00085 there is a signature of Lya absorption sys- 
tems corresponding to a HI cloud at z ss 2.8. 



3.2. Galaxy concentration at z = 2.8 

We noticed that five among 12 galaxies with firm spec- 
troscopic redshifts locate between a narrow redshift range 
from 2.789 to 2.806 (corresponding to a comoving length 
of ~ 17.6 Mpc). In figure [3 we show the spatial distribu- 
tion of 12 galaxies with known redshifts. The five galaxies 
between 2.789 < z < 2.806, shown as filled circles, also 
lie within a ~ 1' X 0/5 sky area (1' corresponds to 1.8 
Mpc (comoving) at z = 2.8). In addition, in the spectrum 
of HDFS00085 at z = 3.16 (shown as a triangle in figure 
EJ) we identified a dump of a continuum around ~ 4600A, 
which coincides with Lya clouds at z 2.80 (see a dip 
around 1 100 A in the spectrum of HDFS0085 in figure^. 
This feature suggests the existence of large amount of neu- 
tral hydrogen at z ~ 2.8 toward the direction. We also 
note that one object among two with probable redshifts, 
HDFS0I295 has an estimated redshift z = 2.798. It is also 
interesting the existence of HDFS00565 and HDFS00526. 
They are a pair of galaxies at z = 2.789, the most bright 
galaxies in our sample galaxies in /-band, and they both 
show a prominent Lya emission, indicating they are expe- 
riencing a massive star formation. Because the comoving 



volume traced by the redshift range 2.789 < z < 2.806 
is less than 2% of the volume surveyed by the range 
2.5 < z < 3.5, the possibility that this concentration is 
caused by chance is less than 10 ~ 4 if we simply assume 
that the distribution of sample galaxies are random. From 
these findings we think that these galaxies would be mem- 
bers of a real galaxy concentration at z ~ 2.8 within a 30 
Mpc 3 or smaller comoving volume. The sky positions of 
galaxies in the redshift range are located at the edge of the 
WFPC2 survey field, and it is possible that this concen- 
tration has an extension out of the surveyed area and the 
number of member galaxies might be larger than currently 
identified. 

Daddi et al. I|2003J1 examined a clustering properties of 
galaxies with photometric redshift estimates in the catalog 
by L03, and found a concentration of galaxies with red 
J — K color (> 1.7) at z p hot ~ 3.0. If we take account the 
possible systematic offset of their photometric redshift at 
z > 2 (see §3.1), the concentration we found coincide well 
with that by Daddi et al. ( 2003 ) , although selected galaxy 
types would be different from each other. 
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Fig. 2. The rest-frame UV spectra of galaxies with redshift identifications. Objects are shown as an order of /si4 
magnitude. The spectra were 3-pixel binned. Units of fx is 10~ 16 erg/s/cm 2 /A. The hatched area represent the 
wavelength ranges dominated by atmospheric lines and where we do not show objects' spectra. The redshifts with 
parentheses indicate that the redshift determination is not definitive. An arrow in the panel of HDFS0085 indicates a 
position of Lya line at z = 2.8, which coincidents with a dump in the continuum (see section 3.2). 



4. Stellar Populations of Star-forming Galaxies at 

z ~ 3 

4.1. Method 

The examination of stellar contents of star-forming galax- 
ies at z ~ 3 has a crucial importance on the understanding 
of their star formation history (SFH), and their connec- 
tion with galaxy populations at different redshifts. Deep 
broad-band photometric data obtained by HST/WFPC2 



(Casertano et al. IMKUl and VLT/ISAAC (L03) covers a 
rest-frame wavelength range from 750A to 5700A for a 
galaxy at z ~ 3 with 7 bands, with which we can con- 
strain some parameters of their stellar contents, such as 
stellar mass, SFH and amount of dust attenuation, with 
reasonable accuracy. In tabled we list the six colors and 
7814 and ifs-band magnitudes for spectroscopic sample 
galaxies. 
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Fig. 3. Same as Fig. |3 for objects fainter than isi4 ^ 24.6. 



We used a population synthesis and evolution code 
PEGASE version 2 (Fioc and Rocca-Volmerange 1997) 
for generating a set of template galaxy spectra. The initial 
mass function adopted is that by Salpeter Q1955[) with a 
mass range of 0.1 Mq - 120 Mq. Two types of SFHs were 
considered, namely, a constant star-formation rate (SFR) 
and an exponentially decaying SFR. 1 PEGASE version 

For exponential SFH, four different values of the decaying 
timescale r (lOMyr, 50Myr, lOOMyr, lGyr) were adopted. In 
order to model different chemical enrichment histories, three 
different SFRs relative to the total amount of baryons within 
a galaxy (consuming 5 x 10 -4 , 5 x 10 -5 , or 5 x 10 -6 of total 
baryons per Myr) were considered in constant SFR models. 



2 takes chemical evolution into account. We started cal- 
culations of synthetic spectra of stellar populations from 
zero-metal gas. The metallicity of interstellar gas of model 
galaxies at fitted ages is sub-solar. Although the measure- 
ments of the metallicity of LBGs still have large uncer- 
tainties, roughly a half Solar metallicity has been sug- 
gested (Pettini et al. I2001f) . Thus we think metallicity of 
our template spectra does not conflict with observations. 
We also calculated spectra of galaxies forming from solar 
initial metallicity, in order to examine the effect of metal- 
licity on derived parameters of stellar populations. Stellar 
mass, ages and dust attenuation obtained using such high 
metallicity models did not change so much from results 
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with low metallicity models; ages were a few dex younger 
and E(B — V) was 0.1-0.2 smaller, but stellar mass did 
not change significantly. So we think metallicity is not a 
dominant source of uncertainty in the estimation of stellar 
populations. 

The procedures for finding the best-fitted models and 
estimating confidence levels we adopted are same as those 
made by previous authors (e.g., Papovich et al. 12001(1 . The 
outputs at 26 time steps from 10 Myr to 2.5 Gyr were 
used for fitting with observed SEDs of 13 sample galaxies 
(11 with firm and 2 with probable redshifts). Before exe- 
cuting the fitting, the absorption by intergalactic neutral 
hydrogen for each redshift of the object was applied, us- 
ing analytic formula in Inoue et al. (2005). We also varied 
the amount of attenuation caused by dust within galaxies, 
from E(B - V) = 0.0 to 1.5 with a step of 0.01, by adopt- 
ing an attenuation law proposed by Calzetti et al. (pOOO). 
For each template SEDs the normalization factor which 
minimize the \ 2 static were explored. The % 2 is defined as 

2 [^ obs fo) Z a X /model fo)] 2 m 

where / D bs(')> /model (i) an d are observed flux, model 
flux and estimated observational error in band i, respec- 
tively, and a is the normalization factor. As observational 
errors we only included photometric errors listed in the 
catalogs. We also tested cases which incorporated addi- 
tional errors possibly associated with systematic differ- 
ences between observing bands. Although the resulting \ 2 
values were reduced in the latter procedure, the changes in 
best-fitting parameters and confidence levels were insignif- 
icant. For objects with Lya emission line, flux densities at- 
tributed to Lya line were estimated from our spectra and 
they were removed from £?45o-band and / or Vesoe-band 
flux densities. For each galaxy we recorded the x 2 values 
calculated for all models, and found the best-fit model 
which returned the minimum \ 2 value. 

The 68% confidence levels were calculated by a Monte 
Carlo resampling. We randomly perturbed the observed 
fluxes, assuming the gaussian distribution of errors. The 
amount of the perturbations were adjusted to match the 
photometric errors. We executed the SED fitting for the 
set of perturbed flux densities, and the resulting SED pa- 
rameters were recorded. The test was repeated 1,000 times 
for each of sample galaxy. Then \ 2 value to contain 68% 
of these simulated data were derived. This value was used 
to determine the range of model parameters for each ob- 
ject with 68% confidence level, by selecting models which 
return % 2 values smaller than that for unperturbed ob- 
ject fluxes. We also tested another way of the confidence 
level estimation with fixed \ 2 values (i.e., without Monte 
Carlo resampling) following Avni 1(1976(1 . The area sizes 
covered by 68% confidence levels are smaller than those 
through Monte Carlo resampling in most cases, while for 
fainter objects the area sizes are comparable. Hereafter we 
use the 68% confidence level areas estimated with Monte 
Carlo resampling as typical errors in our SED fitting. 



4.2. Results 

In Fig. [*"]and[""]we show the results of SED fitting. In ta- 
ble [5J we summarized results in case of constant SFRs and 
initial metal-free gas. Fig. ["""shows the age from the onset 
of the star formation activity and the amount of dust at- 
tenuation (in unit of color excess E(B — V)). The choice of 
SFH (constant or exponentially decaying) does not make 
a large difference in estimated ages and amount of dust 
attenuation for most cases. The median age and color ex- 
cess of sample galaxies are 130 Myr and E(B — V) = 0.29 
for constant SFH, and 80 Myr and E(B - V) = 0.27 for 
exponentially-decaying SFH. When using exponentially- 
decaying SFHs, the best-fit decaying time scales are com- 
parable to the age estimates. The median values of dust 
attenuation we derived are slightly larger than the value in 
some of previous studies {E{B — V) ~ 0.15; e.g., Shapley 
et al. 2001 Stcidcl et al. 1999; but our values are consistent 
with those by Sawicki and Yee I1998*|) . We think the rela- 
tively smaller amount of dust attenuation in LBGs stud- 
ied by Steidel and co-workers are due to their choice of 
filter set and color criteria. Their LBG selection criteria 
in U n — G and G — 1Z color diagram preferentially select 
galaxies with E(B — V) < 0.3 (see figure 2 in Steidel et 
al. 1230*5) 1 . The estimated values of E(B - V) and the dis- 
tribution of our sample galaxies in two-color diagram are 
consistent with each other (see figure ""Q). 

In Fig. the stellar mass estimates are plotted along 
the estimated ages. The difference of stellar mass due to 
the choice of SFH is again small: the median value of stel- 
lar mass is 1.9 x 10 10 Mq and 1.5 x 10 10 Mq for constant 
SFH and exponentially-decaying SFH, respectively. The 
one object which has ~ 10 times smaller stellar mass com- 
pared to other objects is HDFS0271. Because the ifs-band 
magnitude of the objects is about 2 magnitude fainter 
than the average of ifs-band magnitude of other objects, 
smaller stellar mass estimate for this galaxy is reasonable. 

The properties of stellar populations in our sample 
galaxies obtained by the SED fitting analyses are broadly 
consistent with the results of previous studies for LBGs 
at z ~ 3 (e.g., Sawicki and Yee [TUM Shapley et al. I2i""*"""l 
Papovich et al. l2*M)l . 

A critical argument over the estimates of stellar pop- 
ulations in high-redshift galaxies through the SED fitting 
technique is that there might be a missing old stellar pop- 
ulation which is not dominant in the wavelength range 
used for the fitting but contains significant mass of low- 
mass stars. In order to explore such possibility, we also 
tested a SED fitting with model spectra which incorporate 
two components of stellar populations. One is a single-age 
stellar population burstly formed in 2 Gyr ago (z ~ 20 
for an object at z = 3), and the other is a stellar popu- 
lation formed by a constant star formation. The age and 
SFR of the latter were changed in a same way as we did 
in the fitting with a single component. The ratio of total 
baryonic mass (including both stellar and gaseous mass) 
for old and constant components within a model galaxy 
were changed from 1:100 to 10:1. By this test we intended 



I. Iwata et al.: Galaxies at z ~ 3 in the HDF-S 




00 1000 
Age(Myr) 




10 100 1000 

Age(Myr) 



Fig. 6. Distribution of ages and dust attenuation estimates obtained as results of SED fitting for star-forming galaxies 
in the HDF-S. The model spectra of galaxies used for comparison were calculated assuming the constant SFRs (left) 
and exponentially-decaying SFRs (right). Galaxies with probable redshifts (HDFS01295 and HDFS0001) are shown as 
open circles. The error bars show the range of 68% confidence levels. 
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Fig. 7. Ages and stellar mass estimates for star-forming galaxies in the HDF-S. The meanings of symbols and error 
bars are same as Fig. 



to estimate the maximal amount of old stellar component 
which reasonably agrees with observations. For about half 
sample galaxies the best-fitted models were those with less 
than 10% of the total stellar mass was coming from old 
component. For three cases (namely HDFS 0085, 1590 and 
0001) which have stellar mass of the old component larger 
than those from constant SF component, the total stel- 
lar mass becomes 3-5 times larger than the results with 
single component fitting. From this test we confirmed the 
results by Papovich et al. ( 200JJ ; the amount of old stel- 
lar population laying under the young, UV-bright stars 
seems to be not so large in the best-fitted models. Barmby 
et al. H2004J1 reported the first result of SED fitting for 
LBGs using mid-infrared observation executed by IRAC 
on borad the Spitzer Space Telescope, in addition to op- 
tical and near-infrared photometry. They argue that the 
stellar mass and age estimates made with and without 
mid-infrared (rest-frame near-infrared) do not differ sig- 
nificantly. This result also supports the idea that LBGs at 
z ~ 3 have relatively small amount of old stellar pop- 



ulation. Shapley et al. {2005) confirmed the insensitiv- 
ity of best-fit stellar mass estimates to the inclusion of 
rest-frame near-infrared photometry for a sample of UV- 
luminous star-forming galaxies at z ~ 2. However, they ar- 
gued that, by considering the maximal amount of old stel- 
lar components, the upper limit to stellar mass for some 
of star-forming galaxies can be > 5 times larger than the 
best- fitted values with the constant star formation models. 
In our SED fitting results we explored the 68% confidence 
levels of stellar mass constraints in the same way as §4.1. 
The stellar mass of the old stellar components can be 2- 
20 times larger than the best-fitted ones, and in five cases 
old stellar components can dominate > 90% of the total 
stellar mass. Thus although from best-fitted models the 
amount of old stellar component is suggested to be small 
in most cases, our constraints on mass contribution from 
the old component are not strong. 

We should also note that the error in estimating dust 
attenuation is another possible source of uncertainty in 
SED fitting. The best way of estimating dust attenuation 
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would be through the flux ratio of far-infrared emission 
to UV radiation (e.g. Witt & Gordon I3UUU1 Buat et al. 
2005) even for UV-selected samples like LBGs. Burgarella 
et al. I|2005[) performed a SED analysis of nearby galax- 
ies selected by UV or FIR fluxes in Buat et al. 1)2005[) . 
They show that the amount of dust attenuation evalu- 
ated without the far infrared information can be wrong 
by up to 1 magnitude for UV-selected galaxies and up to 
3 magnitudes for FIR-selected galaxies. The uncertainty 
of 1 magnitude on the UV attenuation leads to the uncer- 
tainty of 0.1 magnitude on E{B — V) for the Calzetti's law. 
Moreover, it seems that the dust attenuation law varies in 
the slope and in the strength of the 2175A bump from 
a galaxy to a galaxy. It would enlarge the uncertainty 
of E(B — V). Such a variation of the attenuation law de- 
pending on the galaxy type was examined by Inoue ( 2005) 
theoretically. He suggest that the Calzetti's law is realized 
in a large density medium which is favorable to starburst 
galaxies like LBGs. If it is true, the assumption of the 
Calzetti's law is valid, although we need more investiga- 
tions in this topic, especially with the data of the dust 
infrared emission. 

Shapley et al. : 2001 j examined the correlation between 
the properties of stellar population and spectroscopic fea- 
tures for LBGs at z ~ 3. They found that galaxies with 
young age {t s f, best fit age from the onset of the star for- 
mation, is shorter than or equal to 35 Myr) preferentially 
show broad Lya absorption line and strong low-ionization 
interstellar lines (such as Si II and C n) , while old (t s f > 1 
Gyr) have strong Lya emission lines and weaker interstel- 
lar lines. We examined the relationship between spectro- 
scopic features and parameters of best-fitted stellar popu- 
lation models in our 13 spectroscopic sample. We did not 
find any clear correlation between the Lya emission, Lya 
equivalent widths and the results of SED fitting. The ab- 
sence of any clear trend may be attributed to the small 
number of sample galaxies in our study; the number of 
galaxies with spectroscopic identifications in Shapley et 
al. l|2TifjT|) is 81, which is more than six times larger than 
ours. 

5. Discussion: Comparison with distant red 
galaxies at z ~ 3 

Franx et al. (20(2||) used deep near-infrared images ob- 
tained by FIRES to select galaxies which have photomet- 
ric redshifts larger than 2 and have Js — Ks larger than 2.3 
(Vega) in the HDF-S. They claimed that the red Js — Ks 
colors of these galaxies might be caused by the 4000 A 
break and bright AT-band color would indicate the ex- 
istence of massive stellar population. They also roughly 
estimated the number density of these "distant red galax- 
ies" (hereafter DRGs) and suggested that they might be 
as abundant as half the number density of LBGs. 

We compared the properties of stellar population of 
our sample of star-forming galaxies and DRGs. For this 
purpose we applied a SED fitting procedure for the 13 
DRGs with K < 24.5 in the HDF-S. The method is exactly 



same as that we called "one-component fitting" in section 
4, except that redshifts we used for DRGs were photo- 
metric ones derived by L03. As we noted in section 3.1, 
there might be systematic offsets in photometric redshift 
estimates by L03. To examine the effect of this redshift un- 
certainty in analyses of stellar contents, we executed a test 
in which we adopted redshifts —0.25 smaller than phot-z 
estimates. The estimated amount of dust attenuation did 
not change significantly and for stellar masses, SFRs and 
ages differences were less than 50%, mostly smaller than 
30%. We also tested fitting with "two-components", i.e., a 
maximally old stellar population plus that with constant 
SFRs, and resulting parameters such as stellar masses and 
dust attenuation in young component were not so differ- 
ent from those by one-component fitting. Both constant 
and exponentially-decaying SFHs with various timescales 
were tested. We show estimated ages (from the onset of 
star formation) , dust attenuation and stellar mass in figure 
|H1 As we did for LBG sample, we tested models with solar 
initial gas metallicity as well as those formed from metal- 
free gas, and found that differences in obtained parameters 
were smaller than 68% confidence ranges. As seen in the 
figure, the ages of DRGs are hardly constrained. The 68% 
confidence range of estimated ages for some objects even 
span over an order of magnitude. Such less accurate age es- 
timates for DRGs as compared with those for LBGs would 
be attributed to DRGs' optical faintness. Also DRGs are 
relatively abundant in low-mass stellar population, which 
is less sensible for ages. Given the large uncertainty of 
age estimation, the best-fitted ages of DRGs are on av- 
erage larger than those of LBGs. The median age of our 
DRG sample is 1.6 Gyr, several times longer than that 
of LBGs. The amount of dust attenuation is also system- 
atically larger for DRGs. Although we mentioned before 
(section 4.2), the use of single attenuation law may lead 
the error of ~0.1 mag in the estimation of E(B — V), for 
LBGs best fitted values of E(B — V) are always not larger 
than 0.5, while for DRGs E(B - V) often exceeds 0.5. So 
a systematic difference between E(B — V) for LBGs and 
DRGs seen in figure [3] would be real. 

The errors of estimated stellar mass derived from 
Monte-Carlo resampling were smaller than 0.3 dex, re- 
gardless of large uncertainty of age estimates. It is be- 
cause stellar mass is largely dependent on the fluxes in the 
near-infrared (rest-frame optical) wavelengths. The stel- 
lar mass of DRGs are 10 10 - 10 12 Mq, and the median 
value is 9.4 x 10 10 Af Q for the constant SFR models. DRGs 
have systematically larger stellar masses than our sample 
LBGs. We should beware that, unlike the DRG sample, 
our LBG sample is not complete in A-band magnitude. 
Thus with current data set we cannot exclude the possi- 
ble existence of LBGs which are missed in our sample due 
to the faintness in rest-frame UV and have stellar mass 
comparable to DRGs. On the other hand, SFRs of DRGs 
are fairly large and comparable to those of LBGs: in the 
case of constant SFR models, all best-fit models have a 
few tens to a few hundreds Mq per year. In the case of 
exponentially-decaying SFH, 7 among 13 have on-going 
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SFRs larger than 10 Mq per year and three among the 
remaining six objects have 1 < SFR < 10 in units of Mq 
per year. 

These results of SED fitting imply that many of 
DRGs are not passively evolving galaxies. DRGs would 
be more likely dusty galaxies with relatively large SFRs, 
and they have already accumulated stellar mass more than 
lO lo Af . These results of SED fitting for DRGs are con- 
sistent with some recent studies by other investigators. 
Forster Schreiber et al. (2004) used a sample of DRGs 
from MS 1054-03 field as well as HDF-S for SED fitting 
using template spectra generated with a population syn- 
thesis code by Bruzual and Chariot l|2003|l . The parame- 
ters of stellar populations they obtained agree well with 
ours. 

If the number density of DRGs is as large as a half 
of LBGs' number density, the stellar mass density at 
z ~ 3 based on UV-selected galaxies would be signifi- 
cantly underestimated. Indeed, predictions based on cos- 
mological numerical simulations have suggested the stel- 
lar mass density at redshifts larger than 2 is larger than 
observational estimations using UV-selected sample (e.g., 
Nagamine et al. 12004(1 . However, because near-infrared 
surveys deep enough to detect DRGs are still covers only 
a tiny area of the sky (smaller than several tens of square 
arcmin), uncertainty due to the cosmic variance would be 
large. In addition, they are faint in optical wavelengths 
and spectroscopic follow-up observations are much dif- 
ficult than LBGs. Upcoming wide-field near-infrared in- 
struments would much broaden survey area and increase 
the number of spectroscopically identified DRGs, which 
would be indispensable to derive reliable number density 
of DRGs. 

6. Summary and Conclusions 

We made an optical spectroscopy of star-forming galax- 
ies in the Hubble Deep Field - South with photometric 
redshift estimated to be at z ~ 3. Five firm redshifts 
and two probable redshifts have been newly determined 
in addition to the confirmations of six galaxies with previ- 
ously known redshifts. Six objects among them lie within 
~30 Mpc 3 comoving volume, suggesting the existence of 
a galaxy concentration. 

We investigated stellar populations of sample galaxies 
using optical to near-infrared broad-band imaging data. 
By the comparisons with template spectra we obtained 
best-fit parameters. Typical ages are 20-200 Myr, and 
stellar mass is (0.5-5) x10 10 Mq. These results are consis- 
tent with previous studies on stellar populations of LBGs, 
and we confirmed that stellar content of UV selected star- 
forming galaxies are dominated by young, massive stars 
and dust attenuation is smaller than E(B — V) < 0.5. 
We found no significant correlation between properties of 
stellar populations and UV spectral features, in our 11 
confirmed and 2 possible sample galaxies. 

A comparison of stellar contents in LBGs and DRGs in 
the HDF-S were also made. We found that DRGs tend to 



have larger values than our UV-luminous LBGs for stel- 
lar masses and dust attenuation, and that SFRs in many 
of DRGs and LBGs are comparable. Majority of DRGs 
would be one of the most massive, actively star-forming 
populations at z J> 2. If the number density of DRGs is 
comparable with that of LBGs, current estimates of cos- 
mic star formation rates and stellar mass density at z ~ 3 
based on UV-selected sample would be largely underes- 
timated. To reach a rigid conclusion of the comparison 
between LBGs and DRGs and to draw a complete picture 
of the galaxy evolution in a high- redshift, we need deep 
and wide surveys in near-infrared wavelengths. 
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Fig. 8. Distribution of ages, dust attenuation and stellar mass estimates obtained as results of SED fitting for star- 
forming galaxies (open circles) and distant red galaxies (red filled circles) in the HDF-S. In the left two panels, the 
models assume the constant SFH and in the right panels the models with exponentially-decaying SFH were used for 
comparison with observed SEDs. All calculations of stellar populations in models in these diagrams were started from 
metal-free gas. The error bars indicate the ranges covered by 68% confidence levels. 
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Table 1. Summary of results of spectroscopic observations. The running IDs and designations are taken from the 
version 2 of HST/WFPC2 catalog of the HDF-S (Casertano et al. I2l)l)l)f) . Jgi4 is a corrected isophotal magnitudes in 
F814W filter, also from Casertano et al. I|2UUU|I . and ifs-band magnitudes are 'total' magnitudes (see Labbe et al. 
120031 (L03n. Photometric redshifts (z p h) are from the catalog made by L03. The spectroscopic redshifts obtained by 
our observations are also listed (z sp ). Asterisks denote that the redshifts were determined by our observations for the 
first time. The redshifts with f were based solely on the central wavelengths of Lya lines, and two redshifts embraced 
with parentheses were uncertain due to the poor SN of the spectra. For galaxies with Lya emission lines the rest-frame 
FWHMs and equivalent widths (EW) are also listed (both units are A). 
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00526 


23. 


.293 


22 


.706 


3.657 





.523 


0. 


.185 





.121 





.524 





.007 


2.789 


00059 


23 


.953 


23 


.450 


>2.765 





.937 


0. 


.464 


0. 


.012 





.527 





.263 


2.945* 


00085 


24 


.064 


23 


.192 


1.541 





.658 


0. 


168 





.478 


1 


.134 





.939 


3.170t 


01590 


24 


.156 


23 


.340 


3.211 





.200 


0. 


.112 





.176 





.516 





.039 


2.364 


00153 


24 


.349 


22 


.538 


4.150 





.565 


0. 


192 





.725 


1 


.089 





.181 


2.794 


01295 


24 


.542 


23 


.464 


2.733 





.783 


0. 


.272 





.449 





.813 





.247 


(2.798) 


02376 


24 


.594 


23 


.749 


2.273 





.403 


0. 


163 





.345 





.556 





.077 


2.408* 


00066 


24 


.607 


23 


.076 


2.353 





.656 


0. 


299 





.696 





.912 





.238 


2.801* 


01825 


24 


.643 


23 


.800 


>2.431 


0. 


.830 


0. 


207 


0. 


.250 


0. 


.764 





.525 


3.275* 


00001 


24 


.760 


21 


.916 


1.608 





.444 


0. 


.662 


1 


.722 


1 


.071 





.435 


(2.472)f 


01512 


25 


.228 


22 


.999 


3.402 





.828 


0. 


360 


1 


.161 


1 


.046 





.147 


2.670* 


00271 


25. 


.449 


25 


.432 


1.994 





.253 


0. 


.171 




0.294 


0. 


.248 





.415 


2.806t 



Table 2. Colors of galaxies with spectroscopic redshifts. These colors are derived from magnitude measured with 
0."7 diameter apetures. Isi4 an d Ktot are same as tabled 
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age (Myr) 


E(B - 


V) 


Stellar Mass (M Q ) 


SFR (M Q /year) 


ID 


best 


min 


max 


best 


min 


max 


best 


min 


max 


best 


min 


lilclX 


0565 


50 






0.24 


fl 90 


0.26 


1 7^0-1-10 




9 94.P-HD 
. ^^±en^ J- u 


Q 7Qp i 09 


9 Q4p_i_09 


1 Qp_i_09 


0526 


100 






0.24 


n i q 


0.28 








9 zLOp-I-09 


1 41p_|_09 


^ 79pJ_09 


0059 


50 


1 n 


i nn 




fl 9Q 


0.38 




zl Qzlp-|-0Q 


9 38p_|_1 fl 




1 Qfip_i_09 


fi 9zLp_|_09 


0085 


1400 


200 


2500 


0.02 


0.00 


0.13 


Z: . _■ 1C~ J- \J 


ft Q4p_i_0Q 


4.41e+10 


9 10p-I-01 

Zi . lUCTvl 


1 70p-I-01 


5 5fle+fl1 


1590 


160 


10 


400 


0.20 


0.13 


0.28 


7 45e+09 


1 48e+fl9 


i nfip_i_in 

J_ . UUct -1- U 


t; c;flp_|_fl1 


9 Sflp_i_fl1 


1 QQp_l_fl9 


0153 


500 




i nnn 

1UUU 


n 9Q 


fl 93 


0.36 


a ^9p_i_in 

4:,U^t;T 1U 






1 1 Pip_|_09 




9 '3' : 1p-|-09 


1295 


100 


20 


250 


0.31 


0.24 


0.37 


1 94.p_i_1fl 






1 49p_i_fl9 


QOp+OI 


9 q«pj_09 

Zi . Jub | Uii 


2376 


80 


i n 




98 


fl 1 Q 
u.iy 


0.33 






ft 090-1-00 


ft AOp-UOI 


^ ^^p_|_01 


1 Q^ip_|_09 


0066 


160 


12 


500 


0.32 


0.24 


0.38 


1.95e+10 


4.65e+09 


2.86e+10 


1.44e+02 


6.20e+01 


3.78e+02 


1825 


250 


100 


1600 


0.14 


0.03 


0.19 


7.93e+09 


5.12e+09 


2.10e+10 


3.85e+01 


1.50e+01 


6.70e+01 


0001 


25 


10 


120 


0.52 


0.43 


0.58 


3.33e+10 


1.42e+10 


7.03e+10 


1.39e+03 


5.21c+02 


1.77e+03 


1512 


250 


60 


600 


0.38 


0.30 


0.45 


2.72c+10 


1.54e+10 


3.37e+10 


1.32e+02 


6.15e+01 


3.09e+02 


0271 


80 


10 


800 


0.10 


0.00 


0.17 


8.14e+08 


1.44e+08 


2.47e+09 


1.15e+01 


4.00e+00 


3.15e+01 



Table 3. Results of SED fitting, using template SEDs with constant star formation rates which started with zero 
metallicity gas. "Min" and "max" were lower and upper values for 68% confidence levels obtained through Monte-Carlo 
resamplings. 



